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ABSTRACT 

The Goddard Space Flight Center 2 to 20 keV OSO-8 data on X-ray 
clusters is examined for information which will restrict models for hot 
■ Intracluster gas structures. Our, starting point is the correlations between 
X-ray spectral parameters and optical cluster properties which we presented 
in a previous paper (Mushotzky et al.; Paper I). 

From the correlation between X-ray temperature and velocity dis- 
persion, we conclude that the X-ray core radius should be less than the 
galaxy core radius, if we assume the gas is isothermal. Thus fits of X-ray 
profiles of clusters to isothermal spheres yield a core radius approximately 
equal to the density cutoff radius rather than the true core radius. If 
we instead assume a generalized polytropic structure, from the same 
correlation one may conclude that the gas is more nearly isothermal than 
adiabatic within a few core radii of the center. Therefore we require 
thermal conduction, radiative cooling, or both in cluster atmospheres. In 
both of these interpretations, the Coma cluster and possibly Abell 231 9A 
are quite discrepant. 

From the correlation between bremsstrahlung emission integral and 
temperature, we conclude that hotter clusters contain a larger fraction of 
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their virial mass in the form of gas emitting in the 2 to 20 keV band. 

This conclusion depends on the adoption of the X-ray si-ze measurements of 
clusters from UHURU data by Kellogg and Murray. We find that the "central 
concentration factor" for the emission cannot vary with temperature steeply 
enough to account for the correlation. 

We closely examine relations among the various measures of galaxy 
density and optical cluster morphology. There is evidence for real cor- 
relations between optical central cluster properties and X-ray properties 
which are not "artifacts" of underlying richness correlations as sug- 
gested by Jones and Forman. 

We also predict a population of cool, dim X-ray clusters which 
has not been observed. 

From consideration of the errors in published atomic physics data 
and spectrum calculations, we conclude that the iron abundance deter- 
minations recently quoted for intracluster gas are uncertain by 50 to 
greater than 100 percent from this nonstatistical cause alone. 

Subject headings: galaxies: clusters of— galaxies: -iiitergalactic medium— 


X-ray: spectra 
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I. INTRODUCTION 

In the previous paper (Mushotzky et a1. 1978; hereafter Paper 
I) we presented X-ray results for spatially integrated spectra of twenty 
clusters of galaxies observed by the GSFC Cosmic X-ray Spectroscopy 
experiment aboard the OSO-8 satellite. For these clusters we found 
best fitting temperatures, bremsstrahlung emission integrals, and iron 
abundances (derived from measurements of iron emission line features 
at 6.7 to 6.9 keV). Because our sample of these measured quantities 
is the largest to date, we were able to investigate new correlations, 
between physically interesting X-ray and optically determined cluster 
parameters. In turn, the correlations we found, combined with the 
result that all the clusters could have the same iron abundance, led 
us to confirm that the X-ray emission from clusters in the 2 to 20 
keV band is primarily thermal radiation from hot intracluster gas, 
with perhaps half this gas having come from inside galaxies. In this 
paper we present no further data, but sketch a broader view of the 
significance of our results against the background of other work. 

Our goal is an observational description of X-ray clusters of 
galaxies. Correlations involving the physical quantities that can 
only be determined by spectroscopic measurements contribute much to 
our understanding. In the next few sections we examine certain con- 
clusions that can be drawn from the correlations in Paper I. We find 
from the temperature-velocity dispersion correlation that the X-ray 
core radius should be less than the galaxy core radius, or, alternatively, 
that the polytropic index is about 1.1 for most of our clusters. Analysis 
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of the temperature-emission integral correlation yields evidence that 
more massive clusters accumulate a larger fraction of their mass as 
intracluster gas. We- also reexamine galaxy densities- and optical morphology, 
as they correlate with X-ray properties, for clues as to how mass injection 
by galaxies affects the density structure of the gas. Then, we discuss 
the predictions that can be made from' our correlations and the impact of 
all these results on models for intracluster atmospheres. Finally, we 
devote a section to a critical appraisal of the physical arguments 
used to derive iron abundances from observed equivalent widths -of iron 
line features in X-ray spectra, and we estimate the associated uncertainties 
in abundances derived in this way. Our conclusions are collected into 
§IX. 

It is necessary here to repeat that this discussion assumes that 
the identifications of X-ray sources with the clusters we include are 
correct. Doubt has recently been thrown on *the identification of 
SC1251-28 as an X-ray source (Cooke 1978). In this paper we continue 
our policy in Paper I of treating these twenty identifications as real. 

lie THE VELOCITY DISPERSION VS. TEMPERATURE RELATION 
a) The Observed Relation 

If intra cluster gas is bound in the gravitational potential of 
a cluster of galaxies which are themselves gravitationally bound, the 
virial theorem can be applied to relate the velocity dispersion of 
the galaxies to the cluster mass, and the cluster mass to the X-ray 
temperature. There is evidence that the gas of galaxies in a cluster 
is approximately isothermal, at least in the central regions, in which 
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case the central velocitj^ dispersion Av measures the core mass M 

C* w 

(see Paper I). If some heat transport process also enforces isothermality 

2 ■ 

in the gas, the temperature T should be proportional to Av . If I is 

interpreted as an average X-ray temperature in the emitting region, then 
2 

T o: Av to first order even without isothermality. In Paper I we found 
the best fitting slope of a function t = A Av“ was a = 1,6 + o3, but a 
was pushed below 2 largely due to the influence of the Coma point, 
which is clearly off the trend of the remaining data. We suggested 
a relation leaving out clusters of Rood-Sastry (1971) type cD and B, 
which leaves the other clusters along a line a 2.0, but in fact only 

Coma and A2319A need be omitted from the data set to obtain a ffc 2,0, 

% 

It may or may not be significant that these two clusters are the highest 
temperature X-ray clusters for which T and av have both been measured. 

Coma is the best studied cluster optically, and it seems impossible 
that the central velocity dispersion is significantly higher than 1050 
km s , The velocity dispersion of Abell 2319 is more doubtful. Although 
Faber and Dressier (1976) show some evidence that A2319A and B are 
separate clusters, the X-ray position (Srindlay et 1977) does not 
exclude the possibility that the X-ray emission is associated with 
some combination of A and B. If Av for the combination is adopted, 

A2319 falls on the projected trend of the Av vs. T relation fixed by 

i' 

lower temperature clusters, 

A relation between T and Av very similar to ours was given by 
Mitchell, Ives, and Culhane (1977). Some of their temperature determin- 
ations are significantly lower than ours (notably for Coma and Centaurus), 
for reasons not clear to us. However, even with their lower temperature 
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for Comas Mitchell ^ al_. found that it had a low velocity dispersion 
for its temperature compared to the other cluHers in their sample, 
b) The Extent of the Gas 

Since we have good reason to believe T is roughly proportional 

2 

to Av , we seek an explanation for the value of the proportional ity 
constant. Cavaliere and Fusco-Femiano (1976) have shown that the extent 
of the gas relative to the galaxies is dependent on the ratio of specific 
kinetic energies of galaxies and gas. If m is the average particle 
mass (fully ionized gas with 10% He), this ratio is given by 


Av 


mAv 

kT 


!V 


6.3 


(1000 km s 


-1 


1 keV 
kT 


( 1 ) 


The cluster used by Cavaliere and Fusco-Femiano as an example of the 
value of e is Coma, for which our data give e 0.6 in rough agreement 
with their value. This uses the mean rather than the central velocity 
dispersion. We have seen that Coma lies well below'the trend in the 
Av vs. T relation, and we suggest it may not be a good example. As 
illustration we calculate e and e = (1.17) e individually for each 
cluster in Paper I, Fig. 5, using the most probable value of kT and av 
as given in Tables 2 and 3 of that paper. The results, given in Table 
1, indicate that e ^ 1.1 and K 1.5 for most .clusters. 

When we used the word extent above, we meant the core radius 

1 

of the spherical distribution of galaxies or gas. The core radius is the 
parameter a in the standard approximation to an isothermal sphere, n(r) = 
% (1 + namely a^ for the galaxies and Tor the‘X-ray 

emitting gas. Since the gas density as a function of radius in the 
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outer regions of an intracluster atmosphere depends on the total mass 
of gas, the core radius alone does not determine how far from the center 
gas may be found. The ratio e should be a direct comparison of the 
core radii of galaxy and gas distributions if both are isothermal. 

Values of e > 1 imply that the galaxies are more extended than the 
gas, in terms of core radii. This extent refers to the density distri- 
bution, so the surface brightness, proportional to n through the emission 
measure, will be more sharply peaked, i.e, have somewhat smaller core 
radius. This interpretation of the Av vs. T relation makes it easier 
to understand how intra cluster gas could be largely processed gas, 
ejected from cluster member galaxies and collected near the cluster 
center; if the galaxies were less extended in terms of core radius 
than the gas as has been suggested from previous analyses of X-ray 
emission profiles (Lea 1975, Cavaliere and Fusco-Femiano 1976), some 
heating mechanism for this gas would be required in addition to conversion 
of kinetic energy of galaxy motions into thermal energy. Primordial 
gas falling into the cluster as a continuation of the cluster formation 
process (Gunn and Gott 1972) would also be expected to have a smaller 
core radius than the galaxy core radius, since the gas relaxes much 
faster than the cloud of galaxies. Thus for the combination of roughly 
half processed and half primordial gas implied by our iron abundance 
measurements, we expect e > 1, And as we mentioned in Paper I, the 
best fit to the Av vs. kT relation corresponds to the case where the 
gas is in equilibrium at approximately its gravitational temperature. 
Incidentally, it is difficult to explain the low value of e for Coma 
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by postulating some nonthermal heating of its gas, because radio 
observations Indicate there have been no major outbursts from active 
galaxies in- Coma for about 10 yrs (WrUson 1970)', However, Lea and 

* ' " 'I- 

Holman (1978) show that significant Coulomb heating of the gas can 
occur if the magnetic field is sufficiently low. 

If the distribution of either gas or galaxies is not strictly 
isothermal, then the parameter e has a somewhat more complex interpretation. 
Yet the King (1972) model usually adopted for the galaxy distribution 
is nearly isothermal near the core, and our observations of the strong 
X-ray cluster Perseus indicate that its spectrum is not far from isothermal 
(Smith et^. 1978). 

In any case, values of e > 1 are interesting since fits to X-ray 
profiles of clusters using this standard approximation to isothermal 
spheres in the past (e.g® Kellogg and Murray 1974) have led to X-ray 
core radii a greater than a_. This contradiction is merely the re- 
flection of the fact that the overall gas distribution is not an isothermal 
sphere with infinite extent. If the density (and hence surface brightness) 
falls off rapidly beyond some cutoff radius r , a fit to the surface 

* wU 

brightness using an isothermal sphere will give an apparent a r 

^ X CO 

and not the true a^ which is related to e. Interestingly, the evolutionary 

models of Perrenod (1978a) for intracluster gas would also give apparent 

values of a (in the % to 2 keV band, however) larger than a„, because 

9 

his gas distributions have sharp density cutoffs atr > a , We discuss 

CO g 

these and other models in more detail in §V1I. 
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c) The Polytropic Index of the Gas 

Should one wish to relax the requirement that the gas be isothermal, 

the next simplest equation of state to work with is a polytropic one 

(f « n^ where y is the polytropic index). In spherically symmetric, 

hydrostatic models for cluster atmospheres using this equation, the . 

temperature varies with radius for any 1 < y £ 5/3. For Paper I we 

performed isothermal fits to our spectral data, but we now wish to 

consider models in which the X-ray continuum is the sum or integral 

over contributions from shells at different temperatures. However, 

without fi.tting trial spectra representing these more complicated continua 

(a tedious and statistically unjustifiable procedure)., we can learn 

something about the characteristics of acceptable spectra and hence 

about acceptable gas distributions, although this knowledge will be 

of a rough nature. The technique rests upon the fact that the central 

2 

temperature of hydrostatic models is still proportional to Av , We 
now review enough results on polytropic, hydrostatic, spherically sym- 
metric gas distributions in cluster potentials to enable us to derive 
the equation relating our measured X-ray temperature to av^ and y. 

Defining x as the radius in' units of a , the temperature structure 

y 

is (Lea 1977) 

kT(x) = (x) - ^ (^out^j ^ 

where is the absolute value of the gravitational potential. The 
"outer radius" is determined from 

^ ("out) ™ 


9 


(3) 
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where the "temperature" indicates the degr'ee to which gas is bound 
to the cluster at infinite distance. Thus is the ultimate extent- 
of the gas cloud and the distance at which an exponential drop of the 
radial density law is found. We restrict to negative values for 
this discussion, implying zero gas pressure at infinity and binding 
all gas to the clustero However, this choice excludes cluster gas 
structures connecting smoothly to a universal hot intergalactic medium 
(e.g. Field and Perrenod 1977) or to supercluster gas (Perrenod 1978b). 
The gravitational potential in equations (2) and (3) may be expressed 
with the use of the analytic approximation given by King (1972) for 
the mass distribution, yielding 




iln 


2 hi 3GM 

X + (1 + x^) = “T — . g(x). 


(4) 


where M is the core virial mass and we define the function g(x) for 
future convenience. 

In cases where gas at more than one temperature contributes to 
X~ray emission in the 2-20 keV band, the temperature we have measured 
for each cluster in our sample in Paper I is a particular weighted 
average temperature. This observed temperature parameter, which we 
henceforth call T, characterizes the best fit to the data of isothermal 
trial spectra in the least squares sense. To relate the polytropic 
index y to T, we note that g(x=0) = 1, so that from (2) and (4) the 
central temperature T_ of -our polytropic gas distribution is given 
by 
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1 r 1 

" T{X=0) = ^ I ■ 


(5) 


In order to find how T is related to y, T^, and we performed 

calculations of the expected shapes of continuum spectra that would 

result from polytropic gas distributions and treated them as synthetic 

data sets. We fitted them to isothermal trial spectra- within the same 

energy band (2 to 20 keV) as our real cluster data, using the same 

least squares goodness-of-fit criterion. The results of an extensive 

grid of poly tropic models and isotheraial fits can be expressed as a 

table of values of the ratio T/T^ = t(y, T^, Values of t are 

given in Table 2. It is important to keep in mind that these values 

are dependent on the energy band and goodness-of-fit test. Making 

use of the virial theorem for the core, GM /a = 3(av ) , we have 

^ y w 

numerically 

kT(keV) Z 56 


Av, 


eff 


-1 


1000 km s 


-1 


eff 


(Yeff. T^. Xout) [■'-9(^out>] • 


Given T and Av^ for a cluster, one may find what combination of 
T^, and x^^^ are compatible with the data, by iteratively seeking to 
satisfy equation (6). With an assumption about the unknown value of 
Xout* can obtain limits on Ygff We introduce the notation Ygff» 
the effective polytropic index, to indicate that by this method we 
are estimating the logarithmic slope of whatever equation of state 
P = P(n) is followed by intracluster gas, as long as the gas structure 
meets conditions of spherical symmetry and approximate hydrostasis. 

We discuss this further in §VII. The value of can in principle 
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be determined from spatially resolved observations of clusters, but 
to date there has not been sufficient knowledge of physical conditions 
far from the cluster center to allow it to be determined convincingly. 
Lea (1977) suggests that clusters are observed optically out to x ^ 40, 
implying that one might choose this value, but we have seen that the 
core radius for Isothermal distributions of gas is probably less than 
that of galaxies, and for the gas is nowhere required to be the 
same as the outer radius of the galaxy distribution. Fortunately, 


it turns out 'that when Ygff is determined for some cluster from (6), 
its value is not very dependent on the assumed value of as long 

asx^^^SlO. 

In Figure 1 we repeat the Av and T data from Paper I overlaid 

v» 

with lines of constant Ygff (variable T^) and constant (variable 
Ygjr^) from the results of our synthetic data fits. We emphasize that 

2 

this presentation assumes T <= Av . Estimates of the best value of 
Ygff can be read directly from Figure 1 for each cluster. The arithmetic 
mean of best values for all clusters in the figure is Ygff = 1.08. 

The 90% confidence limits on T restrict the values of Ygff to values 
less than about 1.1 for most clusters, independent of T^, or within 
reasonable bounds. Coma again stands alone as having a higher than 


average Ygff* and due to the crossing of lines of constant Ygff near 
its position in Figure 1, its effective polytropic index is ‘relatively 


indeterminate, 1.15 


Coma 


1 • ■ ■ ■ 

Values of Ygff greater than 5/3 or less than 1 are not unphysical, 
strictly speaking, if thermal instabilities or nonuniform heating occur 
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in cluster atmospheres. We ignore the complications of nonequilibrium 
phenomena and inhomogeneous conditions when we approximate arbitrary 
equations of state by polytropic power laws. 


Abell 2319A lies in a region which is forbidden to equilibrium- models, 
but Y ^ 5/3 is just allowed by the error limits. 

Malina (1978) have derived limits on for Coma and 

Perseus based on available X-ray data at their writing, including spatially 
resolved soft X-ray measuranents. They found 1,18 < y^^^ <1.5 for Coma 
and 1 < <■ 1,15 for Perseus, which agree quite well with our values. 

III. THE GALAXY DENSITY VS, TEMPERATURE RELATION 

In a similar .way we can analyze the correlation we found in Paper I 
between T and the central galaxy density of Bahcall (1977), N^. A 
calibration between and av can be found from Figure 2, where these 
quantities are plotted for clusters having both measurements. Evidence 
that TTq can also be used as a measure of virial mass, a theoretically 
attractive concept, is seen in the appropriateness of the fit 

2 

(galaxies per 4- Mpc^) ^ (22 + 2) ( ?-) ‘ 

° ^ 1000 km s"' 

Using this calibration, the vs, T relation, Figure 7 of Paper I, 

can be analyzed to yield either e or Y^.p.p for each cluster. The numerical 

relation involving Yg.|r.jf is 


kT(keV) J}; 3.7 


^eff “ ^ 
^eff 
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The numbers for most clusters are not very different from what one 
would derive from the Av vs. T relation, but some clusters depart from 
the trend of Figure 2 and thus yield different e and when v$, 

T is used. The results for Coma are significantly different, with 

'xl 1.1 from TTg, substantially lower than Av^. In contrast 

1.2 for A2147 f rom is higher than its value from Av. These 

departures in Figure '2 may be explicable as some combination of optical 
measurement errors in the quantities and Av. Aside from these anomalies, 
the similarity of the best fits to Av vs. T and F vs. T encourages 
the attractive view that .most clusters fall on one joint sequence of 
temperature and virial mass. Measurements of should be inherently 
less sensitive than those of av to errors from inclusion of stray objects 
or exclusion of true cluster members. We' need not emphasize that measure- 
ments of more clusters are needed. 

IVo THE EMISSION INTEGRAL VS. TEMPERATURE RELATION 
The bremsstrahlung luminosity is proportional to the emission 
integral (hereafter El), which is the average of density squared over 

p 

the emission volume, multiplied by that volume, and is denoted <n >V. 

From observations with a field of view larger than one cannot 

reduce El- to its component factors n and V. However, by using all 
available clues to the sizes of the volumes in various clusters, we 
might -hope to make some progress in factoring El. It is- important 
that we try, because the mass of gas cannot be derived from El without 
V and <n>(not <n >). Thus we need some idea of how n depends on radius 


from a -cluster center. 
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The El vs. T relation found in Paper I is particularly important, 
because while T is believed sensitive to the depth of the potential, 

El is sensitive to the amount of gas (if we know the density distribution). 
Thus a correlation between them implies some connection in clusters 
between the mass of gas and the mass of galaxies. The value of such 
a connection is determined by the number of assumptions one must make 
to formulate a useful interpretation. As briefly explained in Paper 
I, if one relates the mass of intracluster gas to the total virial • 
mass of the cluster by (where we must restrict M___ to 

yaS Vii ydo 

represent only that part of the gas which emits X-rays in the 2 to 
20 keV band), then 

El o: , (9) 

2 2 

where f s <n >/<n> > 1 is an overall central concentration factor for 
the gas. The result from Paper I can be expressed as El «cT^, l,6«$p.^ 
5.1, although scatter in the correlation prevents an acceptable fit 
to a single power law. This correlation is probably not due to selection 
effects since only clusters with both low El and low T are strongly 
selected against. The trend can be seen solely from the apparently . 
excluded regions at high EI/lowT and low El/high T. We can eliminate 
the core radius a from (9) right away. Plotting Bahcall's (1975) 
optical measurements of ag for different clusters in our temperature 
sample against T shows no indication that 3g is related to T; in fact, 
she finds a-g ^ .25 Mpc for all clusters with a probable error of +15%, 
an error range smaller than the uncertainties In the T and El data. 


Therefore 
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^ ( 10 ) 

with good confidence. To make further progress we must know how u, 
f, or V depend on temperature. 

The sizes of regions contributing to the 2 to 6 keV flux have - 

been measured by Kellogg and Murray (1974) and all the clusters in 

their list were included in Paper !<> Their sizes are given as apparent 

X-ray core radii a^ for the standard approximation to hydrostatic isothermal 

sphere models’ with infinite extent. Hereafter we call their measured 

parameter R to emphasize that it is an observational quantity and may 
/\ 

not be the true a^ which is related to the density distribution in the 
isothermal region of a realistic model. We treat R as some characteristic 
radius and trust that values of are indicative of relative size 
between clusters. In Figure 3 we show R^vs. T including the upper 
limits in Kellogg and Murray. The four clusters known by their constel- 
lations and having the highest X-ray flux lie on a straight line, whose 
slope is rather dependent on the Virgo point. This line suggests R^^ “ 

1^ with T) 5/4. When Abell 2256 and 262 are included in the fit, 
one obtains instead 1,5 n 1.8, but these two clusters have consid- 
erably weaker flux, and we feel that their sizes are perhaps less certain 
than is indicated by the error limits given in Kellogg and Murray. 

All in all, 1.2 0-6 1.8 is probably a fair representation of the 

slope of the R vs, T relation. 

X 

If we adopt this apparent temperature dependence of R from UHURU 

A 

data (Case A), we have .v^f « T jo subdivide further we consider 
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two simple cases: {Case 1) y fs independent of T, which fixes a range 

of logarithmic slopes for f(T); (Case 2) f is independent of T, which 
fixes the slope of ij(T)o In view of the importance of the Virgo point 
in tying down the low end of the R vs. T correlation and the possibility 

A 

that the gas in Virgo is associated with a single galaxy and has a 
size not comparable with the other clusters, it is also worthwhile 
to consider the consequences of removing this data point. There is 
then no obvious correlation between and T (Case B), and we have 

A 

p R-? 

y f « T from which we can again derive the plausible ranges of 
slopes for f(T) or u(T) in Cases 1 and 2 above. In Table 3 we show 
limits to these slopes. Values of 5 and ? for the four cases have ’ 
little overlap in spite of their wide range. However, for V independent 
of T, 6 or c may be zero, in which event y and f might both be independent 
of T; this null temperature dependence reflects the fact that p may 
be 2 in equation (10). 

In trying to eliminate some of these possibilities, we can say 
nothing a priori about the temperature dependence of y. However, in 
Case lA- the central concentration factor f contains a strong temperature 
dependence. We might ask whether f can vary over the required range 
within the roughly order of magnitude temperature variation of our 
clusters. To be specific, between 2 < T < 13 keV, a factor of six 
in T, we would require f to vary by 6^, 3.2 < 5 < 8,5, or by a factor 
between 310 and 4.1 x 10®. Since polytropic models with y ^ 1«1 
probably good approximations to the X-ray spectral data in Paper I, 
we have calculated values of f for a range of y and and we present 
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the results in Figure 4. The density integrals have been cut off at 
a radius so that f values will be useful for interpretation of 

our EL data. Beyond the temperature is assumed to fall below 

that which contributes bremsstrahilung emission to some given energy 
band. In the 2-20 keV band the lowest temperature that can be detected 
in clusters is about 1 keVo Figure 4 can be used for other bands since 
polytropic models may be normalized arbitrarily. There is a discon- 
tinuity between polytropic models, all of which have the same T(x) 
fall -off with radius for the same x^^^, and isothermal models which 
have no temperature fall -off. Realistic isothermal models have f not 
very different from our y =1.1 calculations. This causes some difficulty 
in the Interpretation of f for low values of y, which can only be resolved 
by considering more physically realistic model So For the present 
discussion, it is sufficient to consider Figure 4. Choosing T^ for 
some model from Figure 1 determines T(xj^g^^)/T^, if TCxj^g^^^) is always 
1 keV. This determines a range, of and from Figure 4a, and 

then a range of f is obtained from Figure 4b. We have calculated f 
for a reasonable range of x^^^^ (from the Abell radius to = 1000) 
and (a range decided on after numerous integrations of models 

giving T agreeing with the Paper I data). Results for a wider range 
of can be deduced from trends in Figure 4. From this procedure 

q c 

we conclude that a dependence of f as steep as r is impossible, 

because f cannot vary that widely in any polytropic model. Even a 
—3 2 

dependence f « T * Is difficult to attain if most clusters have 
y J^.lol. In short. Case lA is improbable for any reasonable variation 
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of Y and from cluster to cluster. 

The case with the fewest objections is Case 2A, since adoption 
of the UHURU size, measurements (keeping in mind the importance of 
Virgo) eliminates Case B, Case 2A implies that hotter clusters have 
' a larger fraction of their virial mass in the form of intergalactic 
gas. emitting in the 2 to 20 keV band. The true temperature dependences 
of V, f, and y may lie between the pure cases we have defined above, 
for example if all three have some dependence, but we suggest that 
some increase in y with increasing T is required. 

The only set of data with which to compare this result is discussed 
by Mitchell, Ives, and Culhane (1977), They consider the effect of 
a restricted energy band on the interpretation of thei r- observed depend- 
ence of L on T, but they Ignore the central concentration factor f. 

A 

r p 

They derive a relation r for Case A and for- Case 

B, which are consistent with our results, if reasonable assumptions 

about f are made. The observational relation on which these M„,^ relations 

gas 

are based is or roughly El <= T^, which lies just within the 

upper end of our range for e. When they plot R from Kellogg and 

o 

Murray (1 974) against their X-ray temperatures, they obtain R « r 

A 

which is slightly too steep to be consistent with our results. 

V. CORRELATIONS INVOLVING GALAXY DENSITIES 

In Paper I we compared Abell's (1958) richness class R with the 
central galaxy density of Bahcall (1977), and we found that cor- 

related better with T and El than R did. These two indicators of galaxy 
density have now been joined by a third, the density <p> calculated 
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by van den Bergh and deRoux (1978) on the basis of R and the cluster 
radii, found by Leir and van den Bergh (1977). Moreover* R, and 
<p> may be related to- the morphological types defined by Rood and Sastry 
(1971; hereafter RS) and by- Bautz and Morgan (1970; hereafter BM). 

In this section we compare all these indicators of galaxy density and 
try to decide which are the most Important for determining X-ray properties 
of clusters. 

First, in Figure 5 we explicitly compare and R for the cluster 
list given by Bahcall (1977), a comparison referred to in Paper I, ' ■ 

The scatter is large, but similar to the scatter in the correlations 
of Paper I that involved N^. In Table 4 we begin to collect a set 
of quantitative measures of scatter, namely mean values and sample 
standard deviations (S,D.) for and <p> among various samples. The 
S,D. of Nq among R = 2 clusters (the most numerous richness class among 
known X-ray clusters) and kT > 6 keV clusters is similar. The latter 
S.D. is larger only because of the influence of Abell 2147 7‘ 

keV, Nq = 12), which is peculiar on other grounds as seen in Figure 
2. Returning to Figure 5, we see that and R are weakly correlated, 
but indeed seem to be measuring different aspects of galaxy density, 
as suggested in Paper I. We have calculated van den Bergh and deRoux *s 
new quantity <p> for the known X-ray clusters in Bahcall 's list. 

In contrast to the previous comparison, the probability of any correlation 
between <p> and IT^ is very small, although <p> should be a more sophisticated 
indicator of galaxy density than plain richness class. This absence 
of correlation can also be seen in Table 4, where the average and S.D. 
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of Nq for two <p>. samples show that <p> fails to distinguish high and 
low ITq, with scatter similar to or larger than in the vs. R comparison. 
Incidentally, <p> correlates equally poorly with T or El for our sample 
of X-ray clusters. Perhaps the method used to determine cluster radii 
by Leir and van den Bergh is too ill-suited for nearby clusters, such 
as most known X-ray clusters. 

Three weak correlations between galaxy densities and optical 
morphological types have recently been suggested, namely between 
and RS type (suggested but undocumented in Bahcall 1977), between R 
and BM type (Leir and van den Bergh 1977), and between <p> and BM type 
(van den Bergh and deRoux 1978). We have given the RS type associated 
with each point in Figure 5, so that the strength of the first of these 
correlations can be examined. More centrally concentrated optical 
appearance does occur more frequently at higher but the scatter 
in Nq among cD and C type clusters (see Table 4) is again large. The 
R vs. BM type correlation is exhibited by Leir and van den Bergh and 
amounts to a small excess of BM type I clusters- in richer classes. 

The <p> vs. BM type relation appears stronger in the table of van den Bergh 
and deRoux, but it is only visible in the mean values of <p> for different 
BM types, since the scatter of <p> within BM types presented in Table 
4 is the largest relative to the mean in the entire table. The redshift 
range considered here, z > .1, is just beyond the redshift of the most 
distant confirmed X-ray clusters and seems to have been chosen because 
it is the range in which the cluster radius measurements of Leir and 
van den Bergh should have the greatest reliability. However, we reluctantly 
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conclude that <p> is not yet a useful quantity for X-ray work. 

For the present purpose, we are left with and R, as well as 
morphological type which seems to have some correlation with N^. But 
we believe that the scatter in the vs. kT, ,vs. El, and vs. 

R correlations is Indicative of a situation in which several factors 
are participating or in which clusters evolve at different rates so 
that we see them in different stages. Let us keep these possibilities 
in mind in what follows. Of the morphological types, it does not seem ' 
to matter much whether RS or BM is used, since the correlation between 
them is fairly good among clusters which have been classified under 
both systems. Since RS type seems correlated with to some extent, 
we concentrate on actual galaxy densities. We have discussed how 
relates to X-ray data rather completely in Paper I, so we now consider 
richness. 

Richness has often been suggested as an important parameter for 
determining the X-ray luminosity of clusters of galaxies. Most recently 
Jones and Forman (1978) and McHardy (1977) have found that increases 
systematically with Abell richness class. Although Jones and Forman 
found.no significant difference between the luminosities of R=0 and 
R=1 clusters, they did find that R=2 clusters were somewhat more luminous 
than the poorer classes. McHardy shows a clear progression of L^‘ through 
R=0, 1 and 2. Perhaps Jones and Forman did not distinguish R=0 and 
1 because the R=0 cluster luminosities were overestimated relative ' 
to R=1 clusters. A comparison of the values of given by McHardy 

A 

and by Jones and Forman, for clusters listed in both of these papers. 
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shows the UHURU luminosities larger than Ariel -5 luminosities by «123 
dex for R=2, by ol05 dex for R=l, but by ,278 dex for R=0 clusters. 

In our own sample, we see no significant increase in El or T with rich- 
ness, but our sample is less than half as large as the samples in the 
above two papers. If we plot T against the logarithm of the galaxy 
counts represented by richness class for our sample (Figure 6) the 
points are scattered' in a roughly circular region around the concentration 
at R=2 and 6 'keV, which appear to be average parameters for presently 
observed X-ray clusters. There is no trend in R vs, T such as the 
trend in vs, T shown in Paper I; the clusters which ruin the pos- 
sibility of such a trend are A1367, A2199, A2319, A2589, and A1254(?), 

In Figure 7 we show the IT^ vs, T correlation restricted to R=2 clusters. 

The correlation persists, showing that and T are correl ated inde- 
pendent of richness. The same result can be seen in the vs, El 
correlation among R=2 clusters. Moreover, if the richness class 0, 

1 and 3 clusters are added to the restricted correlations, in general 

they fall among the R=2 clusters. McHardy (1977) also finds a BM type 

vs. X-ray luminosity correlation which persists when restricted to 

R=2 clusters. However, the converse does not seem to hold: Among 

our cD and B type clusters, no richness dependence of T or El is noticeable. 

The same lack of a richness effect can be seen in Bahcall (1977). 

We are forced to conclude from a critical appraisal of all the available 
data that central cluster properties such as and RS type a're more 
fundamentally correlated with X-ray luminosity (proportional to T^ 

El) than richness is. 
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Jones and Forman (1978) have argued that any correlations 
between morphology and X-ray luminosity are "artifacts'* of a more 
fundamental richness-luminosity .correlation or are selection effects. 

We have shown in Paper I that RS type is better correlated with El than 
R. We have just concluded that (which is related to RS type) is 
more fundamental than.R. Therefore a better view might be that the 
richness-luminosity correlation seen by McHardy and, to a. lesser extent, 
by Jones and Forman .is dependent upon a more fundamental N^-luminosity 
correlation,- The argument of Jones and Forman in support of their view 
boils down to a comparison of luminosity upper limits on cD and B clusters 
to those on R=2 clusters, among those clusters of potentially low luminosity 
because of insufficient observation. It is true that their upper limits 
on the latter are about twice those on the former. Yet we do not feel 
that this situation very strongly favors that existence of low-luminosity 
cD and B clusters oyer the existence of low-luminosity R=2 (or richer) 
clusters. The occurrence of low T, R=2 clusters and low R=2 clusters 
combined with ,our IT^ vs. T correlation indicate that there is more than 
just a richness-luminosity correlation. Perhaps both and R are signifi- 
cant parameters in determining the luminosity, but in different ways. 

Therefore, it behooves us to consider the cases mentioned above, 
with multiple factors or time-varying X-ray properties. One such case 
is that in which central galaxy density, velocity dispersion, and perhaps 
morphology (which reflect the depth and .steepness of the gravitational 
potential) determine the temperature,. while and to some extent rich- 
ness determine the emission integral through a correlation with the 
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total mass of gas. Unfortunately, we have been unable to discover a 
numerical formula involving both and R which predicts the El values 
of our clusters more accurately than the simple power law fitted to 
the Nq vs. El correlation, El R; 4 x 10^^ {Nq/ 20)^ cm“^. If such a dual- 
pafameter formula does not exist, we might look for evidence ,of time- 
varying cluster properties. Courie and Perrenod (1978) and Perrenod 
(1978a) have studied the time dependent response of gas to cluster gravi- 
tational potentials, and they find that cluster atmospheres approach 
structures similar to static solutions fairly quickly and retain these 
structures without large fluctuations after z fb .1 (which is the farthest 
distance to which we see X-ray clusters). However, evolutionary effects 
have been observed or deduced in optical cluster observations (e.g. 

Oemler 1974) and it seems possible that circumstances under which individual 
clusters form could be different enough to allow parameters such as 
El to vary by factors of two or three at the present epoch, in clusters 
which will eventually be closer in El. It is also possible, of course, 
that some of this scatter which worries us is just observational error, 
and we cannot reject that hypothesis until the optical and X-ray obser- 
vations are refined. 

VI. PREDICTIONS FROM CORRELATIONS 
We have not yet used our correlations to make predictions. The 
scatter evident in all of them is a limitation to their utility, and 
all such predictions will have intrinsic uncertainties. Nevertheless 
some correlations are so significant that we cannot resist using the 
best fits at face value. 
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First we convert the known velocity dispersions of clusters into 
a histogram of predicted temperatures using the Av vs. T relation 
(Figure 8a-). A peak in- the distribution of T appears at the surprisingly 
low temperature of 1 to 3 keV, rather than at the average T of clusters 
observed to date (Figure 8c) which is about 6 or 7 keV. We do not know 
what selection effects have prompted optical astronomers to measure 
the velocity dispersions of certain clusters, but at least the number - 
of known av's Is substantial. From a comparison of Figures 8a and 8c 
we predict a population of low temperature clusters, which should be 
of lower luminosity, of lower flux, and harder to observe than hotter 
clusters due to absorption. However, many of these low-velocity-dispersion 
clusters are nearby and should be discovered by HEAO-1 or HEAO-B. 

One can also convert the list of Bahcall (1977) into a similar histogram 
of predicted temperatures from the vs. T relation (Figure 8b). This 
list was compiled from a 2 to 10 keV survey, however, and thus the pre- 
dicted histogram is peaked near the peak of observed temperatures, 6 
or 7 keV, and no clusters below 2 keV are seen. The absence of clusters 
with Nq predicting T > 10 keV is curious and may imply some difference 
between and Av as measures of the gravitational potential for the 
most massive known X-ray clusters. Figure 8c seems to resemble Figure 
8b in its peak and Figure 8a in its breadth. However, the two clusters 
with the highest T in our sample, Abell 1254 and 2142, fall outside 
either prediction of cluster T distributions. In Figure 8c we show 
a Gaussian fit, not including A2142, which has a mean of T 7.0 keV 
and a dispersion a 3.8 keV. With respect to this fit, Abell 1254 
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is 3o from the mean and Abell 2142 is greater than Sa from the mean. 

The isolation of these sources in the temperatl ve histogram may be enough 
to conclude they are not the same type of cluster source as the others, 
perhaps being active galaxies within clusters. 

A prediction of our El vs. T correlation is that if there are 
true cluster X-ray sources more massive than any so far observed, their 
El's will be extremely high, since luminosity rises so steeply with 
temperature. Extrapolation of this relation to the low temperatures 
just predicted above leads to the conclusion that such low luminosity 
objects will be very difficult to find. 

Finally, all our correlations together imply that the observed 
X-ray core radius, obtained from fitting spatial profiles with isothermal 
spheres, is not related to the actual temperature of the gas within 
that radius but is rather some cutoff radius in the surface brightness 
distribution. The reasoning behind this is discussed further in the 
next section. The prediction from it is that observations with better 
resolution than those to date will show a surface brightness distribution 
within the overall profile which implies a core radius smaller than 
the core radius of the galaxies. This prediction apparently does not 
apply to Coma or Abell 231 9A. 

■ VII. IMPACT ON MODELS 

What sort of generalized picture of an Intracluster atmosphere 
is projected by our discussion thus far? The most important aspects 
are; 1) the correlation with depth of the gravitational potential; 

2) the low value of or the high value of e; ^3) the prevalence 



- 28 - 


of iron; 4) ths temperature dependence of u, R , and/or f. From these, 
one visualizes gas coming out of galaxies, falling toward the cluster 
center along with some primordial gas, and building up into a structure 
which is nearly isothermal in its higher temperature regions. Some 
of these points have been made in Paper I, and we refer the reader to 
the final section of that paper. 

Interpretation of the low. values of Ygff requires more discussion. 
Since y = 5/3 represents an adiabatic equation of state and y = 1 an 
isothermal one, we must ask what physical process enforces near isothermality. 
Two possibilities are thermal conduction and central thermal instability. 

Mathews (1978) has made a realistic estimate of thermal conduction 
rates in intracluster atmospheres. The rates depend on the correlation 
length of the intergalactic magnetic field, but if this length is assumed 
to be on the order of the diameter of the galaxies which stir it up, 
the rate is sufficient to enforce isothermality within a few core radii 
of the cluster center. Since conduction across a field is negligible, 
the magnitude of the field Is irrelevant and only its geometry is of 
Interest. Observations of radio tail sources (e.g. Rudnick and Owen 
1977) probably show that the movement of a .galaxy through the. gas causes 
the field to string out behind the galaxy. Turbulence within the tail 
(Pacholczyk and Scott 1976) should have some effect on the correlation 
length long after the galaxy has passed, but whether it should make 
this length much less than 10 kpc is difficult to determine. One process 
which may bear upon the large scale field geometry is the “orbit segre- 
gation” recently suggested by Saslaw (1977). He finds that predominantly 
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radial galaxy orbits accumulate In the center of a slowly contracting 
cluster while predominantly circular orbits spiral in more slowly, due 
to the different time averaged gravitational field sampled by the two 
types of orbits. If galaxies on the radial orbits pull the Intergalactic 
field behind them, the field in the central parts of a cluster could 
also become predominantly radial, leading to an anisotropic correlation 
length. This is further reason to expect a value at least on the order 
of 10 kpc for the radial correlation length. However, observations of 
cluster radio haloes (e.g. Oaffe 1977) may imply that the magnetic cor- 
relation length is small, perhaps only''- 1 kpc. 

The net rate or timescale for heat conduction over some distance 
also depends on the local gas density. Spatially resolved X-ray pro- 
files of clusters prove that the density falls off quickly beyond a 
few core radii (cf. §IIb). It is thus likely that the conduction rate 
drops quickly at about the same radius, causing gas' further out to be 
less coupled to the heat source at the cluster center generated by com- 
pression of the gas and causing the temperature to fall. This is equivalent 
to an increase in beyond this radius. This fact could explain 
why the specific kinetic energy ratio e is greater than 1, indicating 
that the gas is less extended than the galaxies under isothermal assumptions, 
while fits of Isothermal spheres to observed spatial profiles of cluster 
sources seem to show -the gas more extended than the galaxies. The radius 
beyond which conduction is unimportant helps deteinnine the spatial profile 
observed in the kilovolt X-ray range, leading to a discrepancy between 
the actual X-ray core radius and the apparent X-ray core radius for 
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ar) infinite isothermal sphere. The density profiles In Perrenod*s (1978a) 
mass injection models also have cutoffs, but these models are not pertinent 
for other reasons, as discussed below. 

The alternate process capable of causing clusters to have Ygf^ 
close to one is central thermal instability, in which continual central 
collapse causes a steady inflow of intracluster gas. Dissipation in 

i 

the central gas also provides a ready explanation of e > 1, i.e. why 
the gas has a lower remaining specific thermal energy than the specific 

f 

kinetic energy of the galaxies. Realistic models of this type have 
been discussed by Cowie and Binney (1977), Mathews (1978), Cowie and 
Perrenod (1978), and Mathews and Bregman (1978). In this case, we need 
to inquire whether our evaluation of from our synthetic polytropic 
spectra is still valid for these nonhydro static models. Fortunately, 
the gas flow in these models always has a rather low Mach number: This 

implies directly that the advective kinetic energy density in the flow 
is much less than the thermal energy density, so that we do not go far 
wrong in using hydrostatic model spectra compared with data to deduce 
the slope of the equation of state. Hence our conclusion that Ygff 
Jo 1.1 should apply to subsonic dynamic models as well. If one plots 
the pressure and density taken from Figure 1 of Cowie and Binney, the 
logarithmic slope of this effective (radially dependent) equation of 
'"state is very close to one for the range of densities found at 0.5- to 
5 galaxy core radii from the center. Farther out, this equation of 
state steepens, and at the center it flattens, as expected. This particular 
model which they graph does not include thermal conduction, but others 
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in their paper do. It is not hard to. understand why an Inflow model, 
including radiative cooling, produces an isothermal effective equation 
of state without conduction: If gas falls from a radius r.j to rg toward 

the cluster center but cools as it falls, then the gas will not gain 
the full adiabatic pressure increment due to the geometric compression 
factor (r.j/r2) . Instead, since the temperature at r^ is lower than 
the value it would have In an adiabatic flow, the structure of the gas 
can be approximated by a polytrope of index < 5/3, as observed. 

What is the evidence that helps us choose between thermal con- 
duction models, cooling inflow models, or models including both conduction 
and cooling, as an explanation for the low values of Ygff? On the one 

hand, timescales for radiative cooling are not short enough to be of 

• -3 -3 

interest unless central densities are over 10 cm . It is doubtful 

whether densities ever rise that high in all of our clusters with low 
Y and high e. Cowle and Perrenod have indicated that such gas structures 
can be established in realistic situations, but from more sophisticated 
models, Perrenod (1978) seems to reverse this conclusion. Also, the 
fit to the spatial profile of Coma in Cowie and Binney is barely accept- 
able, while polytropic models with a cutoff should provide better fits. 
However, soft X-ray maps of Perseus and Virgo seem to show smaller sizes 
than in the 2 to 10 keV band (Gorenstein et 1977), which may be 

consistent with some cooling at the center. Further spatial observations 

1 

are necessary to resolve this point. On the other hand, thermal conduction 
^In any case, the cooling center would emit largely at lower energies 
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than 2 keV, while the larger region around it emits the 2 to 20 keV 
X-rays observed by OSO-8. This region remains close to a hydrostatic 
polytropic structure with cutof-f, so that cooling -inflow models have 
little effect on the interpretation of the data in Paper I, unless accretion 
onto a central massive galaxy produces a compact source of thermal or 
nonthermal X-rays (but see Mathews and Bregman), 


models require a correlation length for the magnetic field greater than 
about 10 kpc in the radial direction. One cannot show that this require- 
ment is satisfied or not satisfied. Models with both conduction and 
cooling have been discussed by Cowie and Binney, by Cowie and Perrenod, 
and by Perrenod, and should be considered as the most general models 
for analysis. 

The other point requiring more discussion is the interpretation 
of the temperature dependence of El. In §IV we decided that y must 
increase with cluster temperature. Since higher temperature seems to 
be associated with a more massive core (¥ vs. T, Av vs. T),‘ perhaps 

vJ Vi# 

the natural end result of evolution of more massive clusters is higher 
gas densities. Then stripping of interstellar gas into the intergalactic 
medium would be more efficient in hotter clusters, perhaps giving u 
the required temperature dependence. However, it is not clear from 
present theoretical work that a steep enough u(T) would result, and 
our inability to factor y f/V with good confidence forces us to consider 
•other explanations. The key parameter in .understanding the El vs. T 
correlation is the specific mass injection rate per unit cluster mass, 
including galactic and primordial mass injection, since this rate determines 
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how gas builds up around the cluster center. As an example of the com- 
plications possible, note that if a higher injection rate causes f to 
decrease, the effect of greater u is diluted. To make this more clear, 
imagine the evolution of a cluster atmosphere which has already reached 
equilibrium at the center but still is accreting the gas continuing 
to be expelled from cluster galaxies: The central density cannot 

rise above the value it had when the central temperature first reached 
its equilibrium level, unless there is dissipation, and then the rate 
of dissipation determines how fast the density rises. To predict the 
density evolution beyond the core radius and to distinguish the effects 
of continued accretion on the 2 to 20 keV El, it is clear that one needs 
detailed hydrodynamic models with appropriate local physics. 

The territory of detailed hydrodynamic models is still being 
explored. Cowie and Perrenod (1978) and Perrenod (1978a) have considered 
mass injection models with conduction and cooling. Unfortunately they 
do not calculate f or T for their models, so it is impossible to use 
most of their results in our practical analysis. Although Perrenod 
tabulates the X-ray core radius and finds it to be larger than 0.25 
Mpc for all mass injection models except one with saturated conduction, 
he calculates the core radius in the ^ to 2 keV band. The core radius 
in our band should be significantly smaller. The spectra to be expected 
from the nonconductive models of Perrenod may well be inappropriate 
for our clusters both because these spectra deviate from isothermal 
spectra and because the temperatures reached are too high for the 
Coma cluster, which he claims to be modeling. Additional evolutionary 
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calculations such as Perrenod's performed with an eye toward explaining 
our results, should be very rewarding. In particular, the isothermality 
of duster spectra indicates that conduction models should be run with 
a variety of magnetic correlation lengths. 

VIII. RELIABILITY OF THE IRON ABUNDANCE DETERMINATIONS 
Since astronomers are likely to make free use of the result that 
intracluster gas contains about half the solar abundance of iron, it 
is well to consider the reliability of this result. This reliability 
is a separate question from experimental or statistical uncertainty, which 
can be estimated graphically from Figure 2 of Paper I. It is instead a 
question of uncertainty in the underlying chain of physical reasoning 
which links an experimentally determined equivalent width for the 6.7 to 
1 keV iron line feature to an iron abundance. The chain begins with the 
assumptions that the iron line originates in the same gas as the X-ray 
continuum and that this gas is in a statistical equilibrium dominated 
by collisions. The first assumption is wrong if the iron line originates 
in galaxies, but this has been adequately discussed before (references 
in Paper I), or if 'intracluster atmospheres have abundance gradients 
(Fabian and Pringle 1977). The second assumption is wrong if the ionization 
state of the gas is not yet an equilibrium state (very unlikely con- 
sidering the timescale involved for equilibration) or if photoionization 
by a photon field unrelated to the local electron temperature of the 
gas plays a role in determining ionization equilibrium^. If these 

^In models such as polytropic models with large y, the density falls 
off with increasing radius faster than the local X-ray flux from the 
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hot central region, until at some radius, often less than an Abell radius, 
the photoionization rate of elements like silicon is comparable to the 

t 

collisional ionization rate. However, the contribution of photoionized 
gas at such low densities to the overall X-ray flux of a cluster is 
negligible, except perhaps in certain soft X-ray lines. 

assumptions are accepted, the calculation of the iron abundance follows 
from a comparison of the observed llne-to-continuum ratio to the ratio 
expected on the basis of a theoretical calculation of the line emission 
from a hot plasma. For purposes of the following discussion, we will 
call a calculation using atomic data as input and giving an X-ray spectrum 
as output a "spectrum calculation". Uncertainties In these theoretical 
calculations have heretofor been ignored in the literature on X-ray 
clusters, but they have a bearing on the reliability of iron, abundance 
determinations. 

Different experimental groups have derived different iron-to- 
hydrogen ratios from very similar iron line equivalent widths (EW). 

For example, Mitchell ^al^. (1976) have found Fe/H 3 x lO"^ by number' 
from EW ^ 360 eV for the Perseus cluster while we have found Fe/H ^ 

1.4 x 10 ^ from EW ^ 400 eV. The difference is not due simply to different 
measured continuum temperature, but Implies a different interpretation 
of the results of atomic physics available in the literature. Since 
references to atomic physics literature are given in the X-ray cluster 
papers, one could in principle discover the exact reason for any dis- 
crepancy among the scores of individual cross sections used in spectrum 
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calculations. This turns out to be impratical unless each group is 

^ 1 

willing to make the effort to lay all their computer cards on the table. 

It is hoped that this will be done soon, s.ince in. some cases the accuracy 
of the data has surpassed that of the computations. To begin with, 
we will bypass the uncertainties in spectrum calculations relating to 
exactly how the atomic physics .is used and which processes are included, 
arid consider instead the ranges of results within the atomic physics 
itself. 

While line energies assigned to different atomic transitions 
in solar X-ray spectra are probably dependable for astrophysical use, 
excitation cross sections and similar results cannot often be reliably 
computed from solar spectra, because the physical conditions in the 
line formation region on the sun are not well known. An important limitation 
is that density clues are rare and even when found are strictly useful 
only for lines of the one ion employed. In the work used in Paper I 
to derive iron abundances, namely the spectrum calculation of Raymond 
and Smith (1977), use of such solar observational cross sections has 
been almost eliminated for lines contributing to the 6,7 to 7 keV iron 
line blend. In particular, the most important data, collisional excitation 
cross sections for helium-like ions, have been taken from theoretical 
computation rather than solar observations. 

Of course, the theoretical computations of cross sections have 
some range of confidence themselves. Distorted wave approximations 
have been used by a variety of workers, including Jones (1974) and Davis, 
Kepple, and Blaha (1977) and have a potential accuracy of a few percent. 

Since the claimed accuracy of the best calculations is 10 to 20 percent 
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and there Is the possibility of unincluded channels, we favor a pes- 
simistic estimate of uncertainties in such cross sections, roughly 30%, 
This estimate applies to direct colli si onal excitation of helium-like 
transitions, and most other relevant cross sections have larger errors. 
Since the spectrum calculations which use the cross sections also depend 
on an ionization equilibrium calculation, the net uncertainty in EW 
will be somewhat larger, probably 50 to 75%, in the temperature range 
where direct colli si onal excitation is the dominant process contributing 
to the iron line blend. For the cluster iron lines, this range is 3 
key T ^ 30 keV. 

However, line emission for T 3 keV is dominated by satellite 
transitions. Available calculations of the strengths of these satellites 
have greater probable errors than those for lines in helium-like ions. 
Again, comparison of spectrum calculations ostensibly using similar 
atomic input data shows discrepancies symptomatic of the problems; the 
EW vs. T curve given in Figure 2 of Paper I (a representation of the 
results of Raymond and Smith), that given by Bahcall and Sarazin (1978), 
and that estimated by Malina ^ (1978) are similar in shape but 

differ at low temperatures. Until these discrepancies are understood, 
we are left with large uncertainties in iron abundances, a factor of 
2 or more, for lower temperature clusters, for instance Virgo and Abell 
1367. To make matters worse, attention has been drawn to a new process 
by Jacobs et aT_. (1977a, b) which has not been included, to our knowledge, 
in any comprehensive spectrum calculation. This process, autoionization 
leaving an ion in an excited state, impacts the EW vs. T curve at low 
temperatures in several ways.' The changes that must be made in the 
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EW vs. T curve are difficult to predict. 

Finally, it is well to point out that a 10% helium abundance 
Increases bremsstrahlung emission by’ more than 40% over that of a pure 
hydrogen gas, since this emission is proportional to the square of 
the ion charge and since the free-free gaunt factor turns out to be 
larger for helium nuclei than for protons at the same temperature in 
the keV range. Therefore, the choice of helium abundance affects the 
iron abundance derived from the equivalent width* Paper I used He/H - 
8.5% by number. The change in derived iron abundance due to a change 
in helium abundance from 8% to 10% is almost 10%. Of course, the helium 
abundance in intergalactic gas in distant clusters is unknown. The 
inclusion of bremsstrahlung emission due to ions of elements heavier 
than helium is overall a smaller effect. 

In summary, iron abundances for thus far observed clusters of 
average temperature, including all those with T ^ 3. keV, have a probable 
error due atomic physics uncertainty alone of at least 50%, while for 
T -6 3 keV the error increases to at least 100%. 

IX. CONCLUSIONS 

We collect here the conclusions reached from this investigation. 

The first three depend on certain assumptions, which are therefore stated 
again here. 

1) If intracluster gas is hydrostatic or at least has a rather 
low Mach number, and. if the equation of state is restricted to polytropic 
forms, the logarithmic slope of the equation of state is about 1.1 
in the 2 to 20 keV X-ray emitting regions of most clusters'. This con- 
clusion is not very dependent on the level to which gas -fills the potential. 
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as long as ^ 10. 

2) If intracluster gas Is (nearly) isothermal in these regions, 
it has an X-ray core radius less than the core radius of the galaxies. 

3) If the X-ray sizes of clusters measured by Kellogg and Murray 
(1974) and the X-ray temperatures of Paper I are interpreted as showing 
an increase in size with X-ray temperature, then hotter clusters have 

a larger fraction of their virial mass in the form of hot intracluster 
gas than cooler clusters do, 

4) Bahcall’s (1977) central galaxy density correlates with 
X-ray properties of clusters better than richness or the galaxy density 
defined by van den Bergh and deRoux (1978). is related to optical 
cluster morphology. Correlations of and morphology with X-ray properties 
are not "artifacts" of indirect richness correlations. 

5) The probable error in most cluster iron abundance determinations, 
such as those in Paper I, is about + 50% from atomic physics uncertainties 
alone, if the cluster temperature is 3 keV. If T 3 keV, the error 
Increases to more than 100%. 

6) We predict a population of low T, low L clusters that have 
not yet been observed, 

7) The inner regions of intracluster atmospheres must have some 
kind of isothermal! zing dissipation, either thermal conduction, radiative 
cooling, or both. 

8) Therefore the best model for these atmospheres is one in 
which: a) mass is ejected from galaxies; b) it falls toward the cluster 
center; c) it approaches a quasisteady structure which is nearly Isothermal 
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within some boundary radius at which a sharp density- cutoff occurs; 
and d) therefore this cutoff should occur where the observed X-ray profiles 
drop off, while a more sharply peaked profile is entirely within the 
observed profile. 

We would like to acknowledge fruitful discussions with S. M. Lea, 
who suggested calculating Yg-jr-j? f>"om our data, and the useful comments 
of an anonymous referee. R.F.M. is grateful for the support of a 
National Research Council Resident Research Associateship. 



TABLE 1 


VALUES OF THE SPECIFIC KINETIC ENERGY 
RATIO e 


CLUSTER 

• 

AV 

{km-s~^) 

kT 

(keV) 

S 


A401 


1390 

6.7 

1.81 

2.48 

A426 


1396 

6.8 

1.79 

2.45 

.AlO.60 


771 

3.1 

1.20 

1.64 

A1367 


847 

2.8 

1.60 

2.19 

Virgo 


705 

2.2 

1.41 

1.93 

Centaurus 


945 

5.3 

1.06 

1.45 

Coma 


900 

8.9 

0.57 

0.78 

A2029 


1151 

6.2 

1.35 

1.85 

A2U7 


1120 

7.2 

1.09 

1.49 

A2199 


843 

3.2 

1.39 

1.90 

A2256 


1274 

7.0 

1.45 

1.98 

A2319A 


873 

12.5 

.0.38 

0.52 

A2319A,B 


1627 


1.33 

1.82 


. * 


Average of the two values of Av given by Faber and Dressier (1976). 




TABLE 2 


VALUES OF THE TEMPERATURE RATIO 
T E T/T^ 


kT^(keV)' 


3 



10 



30 


^out 

10 

40 

100 

10 

40 

100 

10 

40 

100 

Y 


- 

• 



- 




1.05 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1,00 

1.10 

1.00. 

. .99 

.99 

1.00 

.99 

.99 

.99 

.99 

.98 

1.15 

.99 

.96 

.95 

.98 

.95 

.94 

.97 

.94 

.93 

1.20 

.96 

.92 

.90 

.94 

-89 

.87 

,93 

-87 

.85 

1.3 

.89 

-.81 

.78 

.85 

.75 

.72 

.79 

.70 

.66 

4/3 

.86 

.78 

.75 

.81 

.71 

.68 

.76 

.64 

.60 

1.4 

.82 

.72 

.69 

.75 

.64 

.60 

.69 

.56 

.50 

1.5 

.76 

.65 

.62 

.68 

.56 

.50 

.60 

.46 

.39 

1.6 

.72 

.61 

.57 

.62 

.50 

.44 

-53 

.39 

.32 

5/3 

.69 

.59 

;55 

.60 

.46 

.41 

.50 

.35 

.29 



TABLE 3 


Case A: 
Case B': 


POSSIBLE TEMPERATURE FUNCTIONS 


Case 1 : 

V independent of T, 

■p t’S 

f cc T 


R - r , 3.2 5 ^ 8.5 , 

X 

T1 <v 1 *8 

V independent -0.4 5 ^ 3.1 


Case 2: 

f independent of T, 


1.6^ C ^ 4.3 


“0 .2^ C ^ 1.6 


of T 



TABLE 4 


AVERAGE INDICATORS OF GALAXY DENSITY 


Quantity 

Description 

Number of Clusters 

Average 

SampT 

Averaged 

of Sample 

in Sample 

Value 

S.D. 

^ * 

"o 

R = 2 

14 

26.6 

6.9 


kT > 6 keV 

9 

26.2 

7.8 


kT < 6 keV 

5 

15.2 

3.3 


<p>^ > 60 

8 

22.3 

6.7 

% 

<p> < 60 

8 

' 23.5 

9.0 


El > 8 X 10®^ 

6 

30.8 

3.7 


RS type cD 

12 

25.9 

8.0 

A 

CL 

V 

BM type I, 





0- 1 < z _< 0. 15 

18 

39.8 

29,6 

<p> 

BM type II, 





0.1 < z £ 0.15 

90 

30.7 

18.2 


An samples involving were taken from list of Bahcall (1977). 
Units of Nq are galaxies per ir/4 Mpc^, estimated error is-+ 20%. 

Values of <p> were calculated using expression in van den Bergh and 
deRoux (1978) from data in Leir and van den Bergh (1977). Units of 
<p> are galaxies per unit volume,- normalized arbitrarily in the same 
way as in van den Bergh and deRoux. 
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FIGURE CAPTIONS 


Figure 1 - The central velocity dispersion is plotted against the 

observed best fit single temperature T.' The data, points 
and their errors are repeated from Figure 5 of Paper I and 
can be identified from that Figure. The solid sloping 
lines show the loci of Av^^ and T for polytropic models, 
and the lines are labeled with values of 7 , which is con- 
stant along each line while the central temperature T^ of 
the models varies. The curved dashed lines show the loci 
■ when T^ is constant while y varies. Both sets of lines were 
calculated with = 40. 

Figure 2 - The mean velocity dispersion and the central galaxy density 
(Bahcall 1977) are plotted for all clusters for which both 
quantities have been measured, in order to calibrate 
statistically as an indicator of core mass. The two 
sloping lines cover the approximate range of acceptable 
normalizations for a fit assuming N^ « Av . Velocity 
dispersions are from Faber and Dressier (1976) and Yahil 
and Vidal (1976). Where there is unusual controversy over 
the correct value of av, as for Abell 262 and 2255, two 
values are plotted with the preferred value having a solid 
error line. Abell 2029 is plotted using the average of the 
two values of Av in Faber and Dressier. The typical error 
in quoted by Bahcall is shown in the lower right corner. 

- The X-ray temperature T from Paper I vs, the X-ray core 
radius given by Kellogg and Murray (1974). Errors in T are 


'Figure- '3 



Figure 4 


Figure 5 


Figure 6 


90% confidence except for Abell 262 which is 68% confidence 

and thus plotted with a dashed error line. The sloping dashed line 

is determined by the Virgo, Centaurus, Perseus, and Coma 

points (see text), while the solid line show the flattest 

slope consistent with all. the data .points weighted by their 

errors. Neither line violates any of the upper limits.. 

- The central concentration factor f defined in the text is 

shown in the upper graph as a function of the outer radius 
Xout Several values of y and Values of f have 

been found by integrating the density distributions of poly- 
tropic models out to The lower curve shows the 

temperature at in units of the central temperature T^ for 

several values of The temperature profile in poly tropic 

hydrostatic models is independent of y but depends on x^^^. 

- The list of Bahcall (1977) plotted against Abell richness 

class R, with each point labeled by its Rood-Sastry morphological 
type. There is a correlation between and R at about the 
same level of confidence as between and RS type. The 
scatter in the two correlations is also about equal. 

- Abell richness class R vs. observed X-ray temperature T. The 

right hand scale gives the galaxy counts used by Abell (1958) 
to define his richness classes. Thus the dashed lines determine 
the probable error limits to his galaxy counts, and within the 
dashed boundaries points are plotted at arbitrary levels in 
the vertical direction, their positions chosen to eliminate overlap 



in 90 % confidence error bars on T for different clusters. There 
does not appear to be any correlation between R and T, from 
which we conclude that richness is not a very good measure of 
core mass. 

Figure 7 - The vs. T correlation, shown in Figure 7 of Paper I, when 

only clusters of richness class 2 are plotted. The correlation 
is present at about the same confidence as it is among clusters 
of all richness classes. 

Figure 8 - Histograms of X-ray temperatures: a) predicted from optical 

velocity dispersion lists; b) predicted from Bahcall's 
(1977) IT^ list; and c) best fits from Paper I observations 
by OSO-8. Each temperature T is binned in the 1 keV interval 
determined by its integer part. In a) and b), if T has been 
measured for a cluster and disagrees with that predicted, the 
observed T is binned in a solid numbered box and the predicted 
T is binned in a dashed box containing the same number. The 
clusters having such discrepancies are: 1) Perseus, 

2) A1367, 3) Virgo, .4) Coma, 5) A2029, 6) A2147, 7} A2256, 

8) A2319A, 9) A2199, 10) A2589, 11) Centaurus, and 
12) SC1251-28. In c) a gaussian fit to the observed T 
distribution is superimposed. This fit takes into account 
all clusters except A2142, which was binned as ">20". 
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